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Abstract The excellent spatial resolution of the Chan- 
dra observatory offers the unprecedented possibility to 
measure proper motions at X-ray wavelength with rela- 
tively high accuracy using as reference the background of 
extragalactic or remote galactic X-ray sources. We took 
advantage of this capability to constrain the proper mo- 
tion of RXJ0806.4-4123 and RX J0420.0-5022, two X- 
ray bright and radio quiet isolated neutron stars (INSs) 
discovered by ROSAT and lacking an optical counter- 
part. In this paper, we present results from a prelimi- 
nary analysis from which we derive 2a upper limits of 
76 mas/yr and 138 mas/yr on the proper motions of 
RX J0806.4-4123 and RX J0420.0-5022 respectively We 
use these values together with those of other ROSAT dis- 
covered INSs to constrain the origin, distance and evolu- 
tionary status of this particular group of objects. We find 
that the tangential velocities of radio quiet ROSAT neu- 
tron stars are probably consistent with those of 'normal' 
pulsars. Their distribution on the sky and, for those hav- 
ing accurate proper motion vectors, their possible birth 
places, all point to a local population, probably created 
in the part of the Gould Belt nearest to the earth. 

Keywords stars: neutron . X-rays: individuals 
(RX J0806.4-4123, RX J0420.0-5022) 
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1 Introduction 

The vast majority of isolated neutron stars (INS) are 
being discovered as radio pulsars. However, recent X-ray 
surveys have revealed that a number of these objects 
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manifest themselves at X-ray and 7-ray wavelengths. A 
surprising result of the ROSAT all sky survey has been 
the discovery of a small group of INS exhibiting prop- 
erties at variance from those of the more familiar radio- 
pulsars (see Haberl (2005) for a recent review). First, 
these neutron stars do not exhibit radio emission, al- 
though there are recent claims of the detection at a very 
low radio flux for two of these objects (Malofeev et al. 
2006). Most of them display X-ray pulsations with pe- 
riods close to 10 s. Their X-ray spectra are essentially 
thermal with temperatures in the range of 40 to 100 eV. 
Many ROSAT INSs display shallow broad absorption 
lines on the top of the thermal continuum, which are 
interpreted as proton cyclotron lines or as atomic tran- 
sitions in high magnetic field conditions ( > few lO^'^ 
G). At present, it is still unclear whether or not ROSAT 
INSs constitute a homogeneous group, representative of 
a neutron star population different from standard radio- 
pulsars. The absence of strong radio emission can be due 
either to a very strong magnetic field (in excess of the 
quantum field) or to the fact that the radio pencil beam 
(which narrows at long spin periods) does not sweep the 
earth. Taking into account that only the closest sources 
have been detected by ROSAT because of their intrinsic 
X-ray faintness, their spatial density may well be compa- 
rable to that of normal radio-pulsars. ROSAT INSs may 
be therefore just the tip of the iceberg, a formerly hid- 
den large population of stellar remnants, possibly related 
to the growing number of radio-quiet and X-ray bright 
compact central objects in SNRs (see Pavlov et al. 2004 
for a review). 

The generally low absorption (Nh'^ 10^° cm^^) de- 
rived from the soft X-ray observations suggest relatively 
small distances of the order of a few hundred parsecs. 
Small distances are indeed inferred from the HST par- 
allax of the two X-ray brightest members, RXJ1856.5- 
-3754 (d = 178li7pc, Kaplan 2003) and RXJ0720.4- 
-3125 (d = 330tll°pc, van Kerkwijk & Kaplan 2006). 
Proper motions in the range of 100 to 330masyr~^ are 
now measured for the three X-ray brightest ROSAT INSs 
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(Walter 2001, Motch et al. 2003, Motch et al. 2005). 
These high velocities imply very low Bondi-Hoyle ac- 
cretion rates and thus exclude with high confidence the 
possibility that the neutron star is re-heated by accre- 
tion of matter from the interstellar medium. The ther- 
mal emission of ROSAT INSs is most probably due to 
the progressive cooling of a middle-aged w 10^ — 10® yr 
neutron star. 

In this presentation we report on proper motion mea- 
surements of two ROSAT discovered INSs, RX J0806.4- 
-4123 and RX J0420.0-5022. These two neutron stars 
do not have established optical counterparts which could 
be used to measure their proper motions with optical 
telescopes. Instead we take advantage of the outstanding 
imaging quality of the Chandra observatory to constrain 
their proper motion using two observations obtained in 
2002 and 2005. The first part of the paper describes the 
observations and the method used to analyse them and 
to detect the displacement on the sky. In a second sec- 
tion, we use extensive MARX simulations to validate our 
results. We finally discuss the implication of our findings 
on the origin and distance of the group of thermally X- 
ray emitting but radio-quiet isolated neutron stars dis- 
covered by ROSAT. 



2 X-ray and optical properties of 

RX J0806.4-4123 and RX J0420. 0-5022 

RX J0806.4-4123 was discovered by Haberl et al. (1998) 
in the ROSAT all-sky survey with a PSPC count rate 
of 0.38 cts/s. Follow-up observations with the XMM- 
Newton satellite revealed pulsations at a period of 11.4s 
(Haberl & Zavlin 2002, Haberl et al. 2004). Its thermal 
X-ray spectrum (kT = 96 eV) lies among the hottest 
of this group. As for other ROSAT INSs, the quality 
of the spectral fit is significantly improved by adding 
a shallow absorption line centered around 460 eV. The 
source undergoes very little interstellar absorption (Nh 
= 410^^cm^^ or 1.1 10^" cm~^ depending on the inclu- 
sion or not of the low energy broad absorption line). 
Optical studies have been so far hampered by the low 
galactic latitude of RX J0806.4-4123 {b = -4.98°) and 
no strong constraint exists on the optical brightness of 
the optical counterpart. ESO WFC imaging do not re- 
veal optical counterparts brighter than an estimated R 
magnitude of 22 in the Chandra error circle. Deep radio 
observations failed to detect any source at a position con- 
sistent with that of RXJ0806.4-4123 (Johnston, 2003). 

RX J0420.0-5022 is a slightly fainter INS, also dis- 
covered in the ROSAT all-sky survey (Haberl et al.l999). 
Its PSPC count rate is 0.11 cts/s. XMM-Newton spectra 
show that the source is the coolest of its group with kT = 
45 eV and display evidences of a broad absorption line at 
E ~ 330 eV (Haberl et al. 2004). The interstellar absorp- 
tion toward RX J0420.0-5022 is Nh = l.OlO^^cm-^ or 
2.010^"cm^^, again depending on the inclusion or not 



Table 1 Journal of observations 



Object 




Observation 


exposure 


Roll angle 






date 


time (ks) 


(deg) 


RXJ0806.4- 


-4123 


2002-02-23 


17.7 


322.4 


RXJ0806.4- 


-4123 


2005-02-18 


19.7 


325.6 


RXJ0420.0- 


-5022 


2002-11-11 


19.4 


18.6 


RXJ0420.0- 


-5022 


2005-11-07 


19.7 


23.6 



of a broad line in the spectral model, and is on average 
about twice that of RX J0806.4-4123. Repeated XMM- 
Newton EPIC observations have now confirmed a pulsa- 
tion period of 3.45 s, the shortest of all ROSAT discov- 
ered INSs. A one hour long ESO-VLT exposure suggests 
the presence of a B = 26.6±0.3 mag object in the Chan- 
dra error circle, which owing to its faintness could well 
be the optical counterpart of the X-ray source. 



3 Chandra Observations 

3.1 Detecting proper motions 

We started in Cycle 3 a long term observing program 
aiming at measuring the motion of several ROSAT dis- 
covered INSs, either lacking optical counterparts, or with 
optical counterparts too faint to be measured from the 
ground in a reasonable amount of time. 

A reference astrometric frame can be built from the 
background of AGNs with a relatively high accuracy. 
At low galactic latitudes, a fraction of the field sources 
are identified with remote galactic objects, mostly ac- 
tive stars, which typically have proper motions one or 
two orders of magnitude below that of the neutron star 
and can thus be used to build the reference frame. In 
order to minimize systematic effects, the 2002 and 2005 
observations were acquired with the same instrumental 
setting (aim point, readout mode) and at similar times 
of the year and roll angles (see Table In both cases, 
the neutron star was located at the aim point of the in- 
strument. The duration of the Chandra observation is 
determined by the need to acquire a large enough num- 
ber of well positioned background extragalactic sources. 
In spite of the X-ray softness of RX J0806.4-4123, we 
preferred to use ACIS-I instead of a BI chip of ACIS-S 
to avoid strong pile up and eventual severe degradation 
of the positional accuracy of the INS as a result of the 
distortion of the PSF. For RX J0420.0-5022, which is 
significantly fainter and softer than RX J0806.4— 4123, 
we chose to use ACIS-S to obtain the best position on 
the INS at the expense of a slightly smaller field of view. 

3.2 Source detection 

Data reduction and source detection was performed us- 
ing CIAO 3.3.0.1 and the latest calibration database. The 
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2002 and 2005 observations were reprocessed in an homo- 
geneous way and the best attitude corrections appUed. 
The level- 1 event files were corrected for known process- 
ing offsets and were reprocessed to the level-2 stage. The 
standard Chandra data reduction pipeline randomizes 
the X-ray event positions detected within a given pixel 
in order to remove the "griddcd" appearance of the im- 
ages and to avoid any aliasing effects. In principle, this 
randomization process could slightly degrade the source 
centering acciiracy. We thus had to reprocess the raw 
data in order to remove the pixel randomization follow- 
ing the science thread available at the Chandra X-ray 
center. Source detection was nm on both the normal 
standard processing data and de-randomized data. The 
wavdetect algorithm (wavelet transform) was used with 
pixel scales 1, 2 and 4 well suited to the dc^tecticm of un- 
resolved sources at moderate off-axis distances. We only 
considered sources detected in the central CCDs, chips 
to 3 and 6 and 7 for ACIS-I and ACIS-S respectively. 

The final relative astrometric quality of the reference 
frames also depends on the number of sources common 
to the two observations and on the quality of the deter- 
mination of their positions. We thus tried five different 
source detection thresholds, defined in wavdetect as the 
probability of a false detection at any given pixel. These 
thresholds ranged from 10~* up to a level of 5 10~^ in or- 
der to find the best compromise between the number of 
sources and the mean positional errors, both increasing 
with inci'easing threshold. The largest value implies that 
about 50 false sources per CCD chip enter the wavde- 
tect source list. However, since we only consider sources 
detected in both the 2002 and 2005 observations and lo- 
cated at a maximum distance of one to three arcsec, the 
actual probability that the common source is spurious 
is very low. Best positions are obtained in the energy 
range in which the contrast between the source spectra 
and the unresolved extragalactic and instrumental back- 
grounds is optimal. For the neutron star, measured in 
the same conditions as the reference sources, the choice 
of the energy range is much less critical since the object 
is brighter than any of the background sources. We de- 
cided to test energy bands 0.3-10 keV, 0.3-5 keV, 0.3-3 
keV, 0.5-5 keV and 0.5-2 keV. 



3.3 Matching method 

Our relative astrometric frame determination uses the re- 
constructed equatorial positions of the X-ray sources. We 
allowed for translations of the right ascension and decli- 
nation coordinates as well as for rotation around the aim 
point, resulting from possible errors in the attitude solu- 
tion. A first algorithm simply uses a least square method 
ignoring the source positional errors. In a second step we 
implemented a maximum likelihood method consisting in 
maximizing the quantity; 
L = E,, exp(-i(^)^) 



where dij are the distances between sources i and j 
after transformation and dij the error on their distances, 
both quantities being computed independently for right 
ascension and declination. 



3.4 Simulations 

To our knowledge, the Chandra data base does not con- 
tain repeated ACIS observations of high proper motion 
isolated neutron stars with suitable properties to be used 
as test data. In order to check our detection chain and 
constrain the resulting positional errors, we thus created 
several simulated data sets using the CXC ray-tracer 
MARX 4.2.1 . 

A total of 26 and 12 sources were simulated for the 
fields of RX J0806.4-4123 and RX J0420.0-5022, respec- 
tively. These sources are those common (within 3") to 
the actual 2002 and 2005 observations as detected in 
the 0.5-5.0 keV band and with a threshold of 10~^ for 
RXJ0806.4-4123 and in the 0.5-2.0 keV band with a 
threshold of lO"'^ for RX J0420.0-5022. We assumed a 
power law spectral energy distribution, typical of that of 
the extragalactic population which is known to be domi- 
nant on the X-ray background. Theoretical spectra of an 
absorbed power law of spectral index F = 1.7, undergo- 
ing a hydrogen column density of Nh= 5.21 10^^ cm^^ 
(for the field of RX J0806.4-4123) and Nh= 1.07102° 
cm-2 (for the field of RX J0420.0-5022), were created 
using XSPEC 12.2.1. The column density values were 
derived from the colour excess in the direction of the 
sources as provided by the maps of Schlegel et al. (1998) 
available at the NASA/IPAC Extragalactic Database - 
respectively, E(b_v) = 0.971 and E(b-v) = 0.020 - 
and applying the Predehl & Schmitt (1995) relation be- 
tween the optical extinction and the X-ray absorption, 
Nh= (1.79±0.03)Av 10^1 cm-^. On the other hand, the 
neutron stars were assumed to radiate as blackbodies 
with kT = 96 eV and kT = 44 eV (Haberl 2005). 

Each source was individually simulated giving as in- 
put parameters its equatorial coordinates, photon flux 
and spectral energy distribution, as well as the obser- 
vation exposure and start time (in order to apply the 
proper degree of ACIS hydrocarbon contamination that 
is continuously degrading the quantum efficiency of the 
CCDs, specially at energies below IkeV), the detector 
type at the focal plane and the coordinates of the aim- 
point. The flux normalization of each source was com- 
puted converting the ACIS count rates into flux in units 
of photons s~^ cm~2 in the 0. 1-10.0 keV energy band. 
Events generated for each source were then merged into 
a single event file, which mimics the spatial and flux 
distribution of the actual sources present in the fields 
of RX J0806.4-4123 and RX J0420.0-5022. Background 
X-ray events were extracted from experimental blank 
fields provided by the calibration database. These blank 
fields were reprojected onto the tangent plane of the ob- 
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servations and, for each simulation, the events were ran- 
domly picked from the map so as to reproduce the level 
of background noise of the observation in the 0.5-5.0 keV 
band. Finally, sources and background were merged into 
a single event file. 

The count distribution of the artificial sources is in 
good agreement with a Poisson distribution. We had to 
adjust the input flux of each source simulation to correct 
for the exposure map, vignetting and discrepant counts 
for the few sources located close to CCD edges or bad 
pixels and columns. It also seems that the wavdetect algo- 
rithm systematically underestimates by a few percents or 
more the number of X-ray photons and the input source 
parameters in MARX had to be increased accordingly to 
recover a distribution in source counts similar to that of 
the actual observations. 

A total of 14 pairs of simulated event lists were gen- 
erated for each ACIS I and S detector configurations. 
The main parameters changed were i) the displacement 
apphed to the central source, 0, 0.34, 0.75 and 1", n) 
the offset applied to the 2005 reference coordinate sys- 
tem with respect to the 2002 one, and 0.29" and iii) 
the position of the aim point which was moved by the 
equivalent of 1/4, 1/2 and 2/3 of a pixel in order to 
produce different distributions of the source PSF on the 
CCD pixel grid, although the Lissajous shaped dithering 
should in principle wash out a large part of the effect. In 
each of the 14 configurations, both event lists with and 
without pixel randomization applied were considered. 

We discuss here preliminary results of the ACIS-I 
simulations corresponding to the RX J0806.4— 4123 ob- 
servations. A detailed analysis of all ACIS-I and S simu- 
lations will be presented in a forthcoming paper. 

Owing to its large number of photons, the formal 
errors on the position of the isolated neutron star deliv- 
ered by wavdetect are very small, of the order of 0.02". 
In actual data, however, one should also take into ac- 
count a small additional error which results from un- 
certainties in the transformation between the detector 
pixel and sky coordinates. Deep ACIS-I observation of 
a field in Orion containing over a thousand of bright X- 
ray sources show that this additional error is of the or- 
der of ~ 0.07" (Chandra Proposer's Observatory Guide). 
Therefore, most of the uncertainty on the proper motion 
probably comes from the accuracy with which we can 
link the reference frames of the two epochs. This accu- 
racy is given by = :^:^)J27=iidF? with dP 
being the distance between source i as observed in the 
first and second observation. Most of the configurations 
yield a CTframc close to 0.10". Best results are obtained 
when restricting the energy band in which source detec- 
tion is done to 0.3-5.0 keV or 0.5-5.0 keV. 

In our simulations, the pixel randomization process 
has basically no influence on the accuracy of the frame 
matching and there is no strong dependency either on 
the threshold used for source detection. We show in Fig. 
^the relation between the input displacements and their 




0.2 0.4 0.6 O.B 1 1.2 



Input displacement 

Fig. 1 Measured displacement as a function of the displace- 
ment entered in the simulations. Energy band 0.3-5 keV and 
0.5-5 keV. Filled squares and open squares represent simu- 
lations with and without pixel randomization (small offsets 
applied to the input values for clarity). 

values measured on the simulated data. On the average, 
the motion is recovered with an rms accuracy of the or- 
der of 0.07 to 0.10". In a similar manner, the shifts in 
right ascension and declination are recovered with a rms 
accuracy of 0.07". 



3.5 Results 

As for the MARX 4.2.1 simulations, the (Tframc for real 
data are minimum for the energy ranges 0.3-5.0 keV and 
0.5-5.0 keV for ACIS-I and 0.3-3.0 and 0.5-5.0 keV for 
ACIS-S. The number of sources common to both ob- 
servations varies with detection threshold, energy range 
and maximum allowed distance between the 2002 and 
2005 positions to include the source in the computation 
of the transformation. It ranges from 12 to 39 sources 
for RXJ0806.4-4123 and from 6 to 16 sources for RX- 
J0420. 0—5022. Introducing sources more distant than 1" 
(at large off- axis angles) somewhat worsens the quality 
of the relative astrometry. Detection thresholds smaller 
than 10~^ give results of lower quality as a consequence 
of the decreasing number of real sources detected at the 
two epochs. Using the best combination of energy bands 
and detection thresholds, the mean distance between ref- 
erence sources after transformation is 0.46" for ACIS-I 
and 0.57" for ACIS-S, again with no strong dependency 
on the fact that pixel randomization is applied or not. 
These mean distances are slightly larger by ^ 30% than 
those expected from the formal errors given by wavde- 
tect. A small offset of 0.14" in right ascension and 0.39" 
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in declination is found between the 2002 and 2005 co- 
ordinate systems of the field of RX J0806.4-4123, while 
no really significant translation or rotation is needed be- 
tween the equatorial systems of the two observations of 
RXJ0420. 0-5022. 

For RX J0806.4— 4123, the measured displacement av- 
eraged over the best energy bands and detection thresh- 
olds is ^ 0.1", corresponding to a Icr effect. The error 
on the transformation, crtiamo probably dominates the 
error budget. With a best value of crframc = 0.09" and 
assuming an additional systematic error of 0.07", the 
3(7 upper limit on the total displacement of RX J0806.4- 
—4123 over the 3 years time interval is 0.34" implying 
/i < 114 mas/yr (or /i < 76 mas/yr at the 2a level). 

Because of the smaller central field of ACIS-S and ac- 
cordingly smaller number of reference sources, the mean 
o-framo is 0.19" for the field of RX J0420.0-5022. The 
mean displacement measured is 0.09" corresponding to 
a ~ 0.4cr effect. For RX J0420.0-5022 too, we do not 
measure any statistically significant motion and can put 
a 3(7 upper limit of /i = 207 mas/yr (or /i < 138 mas/yr 
at the 2(7 level). 

Our observing programme has thus established the 
possibility to measure proper motions at X-ray wave- 
lengths with relatively good accuracy. In particular, the 
proper motions of RX J1856. 5-3754, RX J1605.3-I-3249 
and possibly RX J0720.4— 3125 would have been detected 
with ACIS-I and a three year time interval. If Chandra 
lives long enough, very significant constraints or detec- 
tions of the proper motions of the ROSAT discovered 
INSs that are too faint to be observed with optical tele- 
scopes will be possible. 



Table 2 Proper motions, distances and transverse velocities 
of nearby young neutron stars 



Name 




Proper motion 


distance 


Vt 






(mas/yr) 


(pc) 


(kms-i) 


RXJ1856.5- 


3754 


332±1 


178+17 


283 


RXJ0720.4- 


3125 


97±12 


250 


114 


RXJ1605.3-f3249 


144.5±13.2 


< 410 


< 280 


RXJ0806.4- 


4123 


< 76 


240±25 


< 83 


RXJ0420.0- 


5022 


<138 


<340 


< 241 


B0656-H4 




44.1±0.7 


288+^:^ 


60.1 


B1929-hl0 




103.4±0.2 


361+1" 


176 


Vela Pulsar 




58.0±0.1 


287+1? 


79 


Geminga 




170±4 


273±84 


219 




200 400 600 

Transverse speeds (km/s) 



4 Discussion 

The upper limit on the proper motions of RX J0806.4- 
-4123 and RX J0420.0-5022 derived from our Chan- 
dra observations provide significant constraints on the 
space velocities of these two objects. In particular, the 
low value of 76 mas/yr obtained for RX J0806.4-4123 
is well below those measured so far for the three X-ray 
brightest objects but comparable to those of nearby ra- 
dio pulsars (see Table 12). The likely small distance of 
the source {d < 240 pc, Posselt et al. 2006) imply low 
transverse velocities of less than 83 kms~^. 



On Fig. 121 we compare the transverse velocities of 
ROSAT discovered INSs with those of non recycled radio 
pulsars, total and young (age < 3 lO^yr) populations. Ra- 
dio pulsar data are taken from Hobbs et al. (2005). The 
distribution of the inferred 2-D velocities of the ROSAT 
INSs appears consistent with that of the radio pulsars. 
The absence of fast moving ROSAT INSs is intriguing 
but not really statistically significant, since the probabil- 
ity of such a combination remains above the 10% level. 



Fig. 2 Transverse speeds of ROSAT discovered INSs com- 
pared to those of non recycled radio pulsars. Data taken from 
Hobbs et al. (2005). The black thick line and the thin ma- 
genta line show the histogrammes of the entire population 
and of the young (age < 3 lO^yr) respectively. Measured val- 
ues and upper limits for the ROSAT INSs are shown on the 
top. 



According to Hobbs et al. (2005), young radio pulsars 
have a mean 2-D velocity slightly larger than old ones. 
However, the two distributions overlap, as seen on Fig. [3 
and the significance of this difference remains uncertain. 
The ATNF catalogue (Manchester et al. 2005) contains a 
total of 6 radio pulsars with d < 1 kpc, younger than 4.25 
Myr and with measured proper motion. They have 2-D 
velocity in the range of 60 to 354 kms-i with a mean of 
177 kms~^, clearly consistent with that observed for the 
ROSAT INSs. 

The measurement of proper motions can potentially 
provide a wealth of important informations on the evolu- 
tionary status, origin and age of the neutron star. A large 
spatial velocity, typically in excess of 30 kms^^ rules out 
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the possibility that the star is reheated by accretion of 
interstellar material. 

In this respect, the upper limits on the 2-D velocities 
of RX J0806.4-4123 and RX J0420.0-5022 are still large 
and the possibility that they accrete from the interstel- 
lar medium cannot yet be seriously considered. The fact 
that these two INSs share many properties such as X- 
ray energy distribution and pulsation periods with high 
spatial velocity ROSAT INSs also argues for a similar 
nature. 

If measured with enough accuracy, the proper motion 
vector of the neutron star can also be used to compute 
possible past trajectories assuming a range of present dis- 
tances and radial velocities. OB stars progenitors of neu- 
tron stars are not distributed randomly on the sky but 
rather concentrate in relatively large associations shar- 
ing a common space motion. The exact fraction of early 
type stars found in OB associations remains in many 
cases uncertain due to the errors on the photometric dis- 
tances and the generally fuzzy boundaries attributed to 
the associations. Garmany (1994) shows that the relative 
fraction of O and B stars found in associations is larger 
than 75% and 58% respectively. The fact that a group 
of backwards trajectories crosses the past position of a 
known OB association is thus a hint, no definite proof 
since other trajectories are possible, that the neutron 
star progenitor was located there when it underwent su- 
pernova explosion. It becomes then possible to compute 
possible travel times which can be compared to theo- 
retical cooling curves, with the spin down time and, if 
associated to a SNR, with the expansion time. For seven 
nearby and young isolated neutron stars listed in Table 
13 a tentative birth place has been proposed in the recent 
literature. 

We show on Fig. 13 the positions projected on the 
galactic plane of neutron stars younger than 4.25 Myr 
and located at distances less than 1 kpc. Radio pulsar 
data were extracted from the pulsar catalogue main- 
tained by the Australian Telescope National Facility (Ma- 
nchester et al. 2005). We also show the boundaries of the 
classical OB associations as defined in Humphreys (1978) 
and de Zeeuw et al. (1999). For the Gould Belt, we use 
the determination of Perrot & Grenier (2003) based on 
HI and H2 clouds and Hipparcos distances to the nearby 
OB associations. 

Neutron star travel times from their supposed birth 
place to their present location are difficult to determine 
with great accuracy. The large extents of the possible 
parent association and the unknown radial velocities usu- 
ally yield rather large uncertainties. In the case of the two 
youngest nearby radio pulsars, the Vela pulsar (lO^'yr) 
and B0656-I-14 (lO^yr), there is a good agreement be- 
tween the estimated travel time, the pulsar age com- 
puted assuming magnetic dipole braking and the age 
of the remnant of the parent supernova (Hoogerwerf et 
al. 2001; Thorsett et al. 2003). The absence of obvi- 
ous supernova remnants at the possible birth place of 



B19294'10 and Geminga is consistent with their older 
ages. In this respect, the lack of conspicuous remnants 
in the surroundings of the hypothetical birth locations of 
the ROSAT discovered INSs is also consistent with the 
relatively large ages (> 3 lO^yr) derived from the flight 
times (see Table OJ. 

According to Perrot & Grenier (2003), the Gould Belt 
has an elliptical shape with major and minor semi-axes 
of 373 ± 5 pc and 233 ± 5 pc respectively. The Sun lies in- 
side the Belt and is located about 104±4pc offset from 
its center. It is currently closer to the region of the Gould 
Beh populated by the Sco 0B2 and Vela 0B2 + Trum- 
pler 10 associations (see Fig.O). Backwards trajectories 
suggest an origin in the closest region of the Belt for 
RXJ1856.5-3754, RX J0720.4-3125, RX J1605.3-f3249 
and B1929-I-10. Two INSs, Geminga and B0656-I-14 were 
possibly born in or close to the Orion complex, on the 
opposite side of the Gould Belt. As visible in Fig. |3| all 
ROSAT discovered INSs lie in the half sky centered on 
Sco 0B2 and containing the closest part of the Gould 
Beh. 

No such sensitivity anisotropy exists in the ROSAT 
all sky survey at the level of 0.14 cts/s, the PSPC count 
rate of RX J0420.0-5022 which is the faintest of the 
ROSAT discovered INSs (Voges et al. 1999). In order 
to detect and recognise the extreme softness of such a 
source a minimum of ^ 200 s exposure time is needed. 
Only 10% of the whole sky had exposures shorter than 
this value. Furthermore, the distribution of the low ex- 
posed regions are not specifically directed towards the 
far side of the Gould Belt where ROSAT INSs are ap- 
parently not detected. 

The very soft thermal-like energy distributions of the 
ROSAT INSs render their detectability very sensitive to 
the amount of interstellar absorption on the line of sight. 
The distribution of the local interstellar medium has 
been the subject of extensive studies, mostly based on 
the measurement of the NalD doublet at 5890A (Welsh 
et al. 1994; Sfeir et al. 1999). Our current environment 
is rather unusual since we appear to be located in a local 
bubble cavity of very low interstellar gas density. The 
bubble is surrounded by a "wall" of interstellar matter. 
Its shape is very irregular with the closest distance to 
the edge of the bubble being at ~ 60 pc in the direction 
of the Galactic center. Being based on more than a thou- 
sand of directions, the Nal absorption map of Lallement 
et al. (2003) provides a good description of the local in- 
terstellar medium up to distances of about 350 pc for 
absorptions up to a few 10^° cm~^. Their map shows 
no evidence for any asymmetry which could explain the 
absence of detection of INSs in the half sky where they 
apparently miss. 

Therefore, although based on small number statis- 
tics, the distribution of ROSAT discovered INSs on the 
sky might well reveal the major role played by the lo- 
cal Gould Belt in their generation. Popov et al. (2003) 
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Fig. 3 Positions of nearby neutron stars and OB associations projected on the Galactic plane. Red boxes show the OB asso- 
ciations boundaries. Hipparcos stars with a probability higher than 75% to be linked to the OB associations are represented 
by black dots. Blue lines and filled squares show the possible positions of the ROSAT discovered INS, assuming a distance 
range of 100 to 400 pc for those which do not have distance estimates. Radio or 7-ray pulsars younger than 4.25 Myr and 
located within 1 kpc are shown as magenta symbols (when a parallax distance exists) or green symbols when the distance is 
estimated from dispersion measurements. The Gould Belt is shown as a thick dashed line. 
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Table 3 Possible birth places of nearby young neutron stars 



Name 


spin down Age 


Travel time 


Birth place 


Reference 




(yr) 


(yr) 






RXJ1856.5-3754 


7 


4 10" - 1 10" 


Upper Sco OB2 


1,2 


RXJ0720.4-3125 


210^ 


6 10" - 3 10" 


Tr 10 + Vela OB2 or Lower Sco OB2 


3,11 


RXJ1605.3+3249 


? 


~10" 


upper Sco 0B2 


4 


B0656+14 




10'^ 


Monogem Ring 


5,6 


B1929+10 


310'' 


1 10" 


Upper Sco OB2 


7,8 


Vela Pulsar 


10" 


10" 


Vela OB2 


9,8 


Geminga 


3.410^ 


3.4 10" 


Cas-Tau or Ori OBI 


10 



1) Kaplan 2003; 2) Waher 2001; 3) Motch et al. 2003; 4) Motch et al. 2005; 5) Brisken et al. 2003; 6) Thorsett et al. 2003; 7) 
Chatterjee et al. 2004; 8)Hoogerwerf et al. 2001; 9) Dodson et al. 2003; 10) Pellizza et al. 2005; 11) Kaplan & van Kerkwijk 
2005. 



have indeed highlighted the importance of the Gould 
Belt for understanding the bright part of the observed 
X-ray logN(>S)-log S INS curve. 



5 Conclusions 



The high astrometric quality of the Chandra Observa- 
tory offers the possibility to measure proper motions of 
X-ray sources with an unprecedented accuracy. Our ob- 
servations and simulations show that displacements as 
small as ~ 0.3" can be detected. Applying this method 
to two ROSAT discovered INSs, we derive an upper limit 
on the transverse speed of RX J0806.4— 4123 at the lower 
end of the observed distribution for radio pulsars while 
less constraining limits are obtained for RX J0420. 0—5022 
We show that the spatial distribution and possible birth 
places of this particular group of INSs argue in favour 
of a local population and for a production dominated by 
the closest part of the Gould Belt. 

Acknowledgements We acknowledge the use of the ATNF 
Pulsar Catalogue available at ,http:/ /www.atnf.csiro.au/| re- 
search /pulsar /psrcat . 
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